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Low temperatures are necessary for the observation of strongly correlated quantum phases of fermionic atoms 
in optical lattices. We analyze how the temperature of a Fermi gas is altered when the fermions are loaded into an 
optical lattice with an underlying harmonic confining potential and show how the temperature can be measured. 
The temperature of the atoms in the optical lattice determines the fraction of doubly occupied lattice sites of a 
two-component Fermi gas. We analytically calculate this quantity and find a strong temperature dependence. 
This fraction can be measured by studying the production of molecules in the lattice using a Feshbach resonance 
which allows for precise thermometry of atoms in an optical lattice. 
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Ultracold fermionic atoms in optical lattices are a promis- 
ing tool to simulate problems from condensed matter physics. 
Interesting phase transitions between many-body quantum 
phases in the lattice are predicted JH0I- F° r these phase tran- 
sitions the temperature is one of the key control parameters 
and a precise knowledge of the temperature of the atoms in 
the lattice is necessary to experimentally determine the phase 
diagram. Proposals for cooling the atoms in the lattice in or- 
der to reach interesting quantum phases such as the BEC-BCS 
transition in lattices or an anti-ferromagnetic phase have been 
devised (jHHEI- However, standard methods for thermome- 
try - such as observing the rounding-off of the Fermi edge [5] 
- fail in the inhomogeneous lattice since they turned out to be 
dominated by the trapping potential rather than by tempera- 
ture 00. 

Previous theoretical investigations of ultracold fermions in 
optical lattices have often neglected the role of the harmonic 
confining potential. For ongoing experiments, this confining 
potential is an unavoidable ingredient and its consequences, 
such as an inhomogeneous filling of the lattice and the co- 
existence of insulating and conducting regions, must be con- 
sidered. In a typical experimental sequence a quantum degen- 
erate Fermi gas is initially prepared in a purely harmonic opti- 
cal dipole potential and subsequently transferred into a three- 
dimensional optical lattice with underlying confining potential 

0. 

The transfer of a Fermi gas from a harmonic trapping poten- 
tial into an optical lattice leads to a significant modification of 
the density of states. A change in the density of states can re- 
sult in both an adiabatic cooling or an adiabatic heating of the 
gas. We find that for the situation realized in the experiments 
the gas is adiabatically heated which poses a difficulty for 
future experiments aiming at very low temperatures. To de- 
termine the temperature in the lattice we analyze the fraction 
of doubly occupied sites of a noninteracting two-component 
Fermi gas 1 8]. This can be measured experimentally by study- 
ing the formation of dimers in the optical lattice on a Feshbach 
resonance |9] or using photo association We cal- 

culate the fraction of doubly occupied sites analytically and 
find a strong temperature dependence which makes it ideal 
for thermometry of ultracold fermions in an optical lattice. 

A noninteracting two-component Fermi gas in an optical 
lattice with an underlying harmonic confining potential can 



be described by the Hamiltonian 

<j,l> 3 

Here a — {f , J.} denotes the spin of the particle, t denotes the 
hopping matrix element between adjacent lattice sites, Cj a and 
Cj a are the creation and the annihilation operators on the lat- 
tice site j, respectively, and nj a — c^Cja- is the number oper- 
ator. Furthermore, m denotes the atomic mass, uo the trapping 
frequency of the harmonic confining potential and d = A/2 
the lattice spacing given by the laser wavelength A. The in- 
teractions between atoms in different spin states are assumed 
to be zero, as realized experimentally by using a Feshbach 
resonance 0]. The trapping frequency u> of the confining po- 
tential is determined by the potential depth V and the waist 
w of the Gaussian intensity profile of the laser beam. For a 
three-dimensional lattice generated by three mutually orthog- 
onal laser beams of identical waist w and potential depth V 
one obtains ^uj 2 = ^C. 
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Let us first consider the single particle density of states 1 13] 
in the low-tunnelling limit in a three-dimensional optical lat- 
tice. There the particles are localized in the potential wells and 
we restrict ourselves to a single band approximation which is 
justified by the experiment Q]. For determining the density 
of states we assume the width of the lowest Bloch band At to 
be smaller than the energy offset between neighboring lattice 
sites due to the harmonic confinement. This approximation is 
valid except for a very small region at the center of the har- 
monically trapped cloud, where delocalized atoms should be 
considered. However, only few atoms are affected and we 
neglect their contribution in the analytical calculation. The 
presence of the harmonic confinement modifies the density of 
states significantly: for a particle on the lattice site j, the po- 
tential energy is ^to 2 (-|) j 2 . Allowing for only one particle 
per lattice site, the density of states in three dimensions reads 

P3D (E) = 2-kE 1 ' 2 (s^y V2 . (2) 

This shows that the density of states of fermions localized 
in an optical lattice with underlying harmonic confinement 
has the same power law dependence as that of free fermions. 
This analogy can be understood by the following argument: 
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in free space the fermions are characterized by their momen- 
tum k and the eigenstates are equally spaced in momentum 
space. Each fermion has a kinetic energy E oc k 2 . In the 
low-tunnelling limit with harmonic confinement the fermions 
occupy equally spaced lattice sites at positions jX/2 and have 
a potential energy E oc j 2 (A/2) 2 only, which results in the 
same dependence of the density of states on E as in the free 
fermion case. From the density of states we directly obtain the 
Fermi energy in the tight-binding limit: 
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In one and two dimensions we find the same analogy to 
the free fermion case, namely P2d{E) = 8tt / (mu> 2 A 2 ) and 
p XD {E) = (mLu 2 X 2 E/8)-^ 2 , or generally p(E) oc E d l 2 ~ x 
with d being the dimension. 

To estimate the validity of the low-tunnelling approxima- 
tion we compare the analytical results with a numerical calcu- 
lation where we diagonalize the Hamiltonian Q exactly. We 
compute the density of states and fit the result with a function 
p(E) — a ■ E v to obtain the exponent v. For a large poten- 
tial depth, the value of v approaches 1/2, as expected from 
equation (0. Similarly, for a two-dimensional quantum gas in 
an optical lattice the exponent approaches 0. The calculations 
were performed on 300 2 lattice sites in 2D and 100 3 lattice 
sites in 3D. 

When atoms are loaded into an optical lattice the density 
of states changes from the three-dimensional harmonic oscil- 
lator result p(E) oc E 2 to the low-tunneling results p(E) oc 
Adiabatically changing the density of states of a sys- 
tem can be used to cool or heat a gas, provided that collisions 
keep the gas in thermal equilibrium 1141 1 1 511 . We assume an 
adiabatic evolution of the gas during the loading procedure 
which implies that no heat is exchanged with the environment. 
The interactions between the atoms are supposed to be suffi- 
ciently large to ensure thermalization of the sample but small 
enough so that the energy spectrum of the gas is not altered. 
In such a situation the gas evolves along a path of constant 
entropy S. The entropy S of a Fermi gas is given by flcjll 
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with the total energy £ = Ep(E)f(E)dE, the chem- 
ical potential p and the temperature T, n labels the en- 
ergy eigenstates E n and f(E) = (e < - E ~^/ kBT + l) 1 is the 
Fermi distribution function. The chemical potential is de- 
termined from the normalization to the total particle number 
N = p(E)f(E)dE. For low temperatures the entropy 
in the two limiting cases can be obtained using the Sommer- 
feld approximation. For noninteracting fermions trapped in 
a three-dimensional power-law potential V(r) oc r m 



the entropy is given by 1117111 

— = Ntt 2 - + - 
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We now compare the result of the entropy for a given number 
of particles in a harmonic trap (a = 2) with an ideal Fermi gas 
in an optical lattice with underlying harmonic confinement. 
As discussed above, the latter has the same density of states 
as free fermions in a box potential and therefore corresponds 
to a = oo. Consequently, the entropy in a harmonic trap is a 
factor of 2 larger than for the fermions in the optical lattice. 
Thus when the transfer from the harmonic potential into the 
lattice is performed adiabatically it is accompanied by a tem- 
perature increase by a factor 2. The effect of the turn-on of 
the lattice on the atom number statistics at the lattice sites has 
been calculated for a one-dimensional Fermi gas fl^l . 
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FIG. 1: Curves of constant entropy in the lattice. The calculation is 
performed by diagonalizing Hamiltonian Q with N = 50000 atoms 
trapped in three-dimensional optical lattice with the waist of the lat- 
tice laser being w = 70 /im, according to the situation in |9||. The 
values shown at zero lattice depth are the harmonic oscillator result 
without lattice. For s = V/E r < 5 the tight-binding Hamiltonian 
may give incorrect results and therefore we refrain from numerical 
calculations in this regime. 

We numerically calculate S for various temperatures and 
lattice depths (see fig. 0, similar to what has been done pre- 
viously for the homogeneous lattice |;3?]. For a large potential 
depth the numerical calculations agree with the analytical pre- 
diction of equation l|5}- Moreover, when moving along a line 
of constant entropy as the lattice depth is increased we observe 
that the temperature never drops below its initial value in the 
harmonic oscillator trap. This implies that no cooling can be 
achieved during the transfer into the lattice. 

The quest for obtaining ultra-low temperatures in the op- 
tical lattice is accompanied with the need for precise ther- 
mometry. We show that the temperature can be obtained from 
the fraction of doubly occupied sites of a noninteracting two- 
component Fermi gas. A similar problem has been studied 
numerically for bosonic 1 20] andfermionic |01 atoms. We an- 
alytically calculate the number of doubly occupied sites N2 
assuming that both species are distributed according to the 
same Fermi distribution function f(E) = f]{E) = fy(E) 



(5) 



N 2 = 



p(E)f 2 (E)dE 



(6) 
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In two dimensions (i.e. for v = 0) simple analytic expres- 
sions are obtained. We consider a Fermi gas strongly confined 
along one axis to create a two-dimensional quantum gas which 
is subject to two crossed optical lattices in the transverse di- 
rections. The chemical potential as a function of temperature 
is fi = k B T In (e EF t kBT — l) and the fraction of doubly oc- 
cupied sites is 




Temperature T7T F Temperature T/T F 



For a density of states p oc E v the fraction of doubly occupied 
sites is given to first order in fc^T / Ep by 




This expression shows that there is a linear dependence of the 
fraction of doubly occupied sites on temperature in the degen- 
erate regime k B T -C Ep. The slope depends only on the 
exponent of the density of states, but not on lattice depth (see 
figure^. This makes this quantity ideally suited for thermom- 
etry. 

For the special case of a deep three-dimensional optical lat- 
tice with v — 1/2 (see equation (0) we use the known ex- 
pression for the chemical potential n = Ep ^1 — -^(^eT") 2 ) 
and expand n-x to third order in k B T/Ep\ 




FIG. 2: a) Comparison of the analytical formula {7) (solid line) with 
the numerically calculated double occupation of a two-dimensional 
quantum gas in an optical lattice for various parameters by diago- 
nalizing Hamiltonian Q. Squares: s=25, w=70 pm, iV=5000, stars: 
s =15, w=100/xm, iV=3500, triangles: s=15, w=50 pm, N=2500. 
b) Comparison of the analytical formula il 21 (solid line) with the 
numerically calculated occupation in a 3D optical lattice for vari- 
ous experimentally relevant parameters by diagonalizing Hamilto- 
nian G). Circles: s=15, w=50 pm, iV=25000, triangles: s=20, 
w=lQpm, iV=15000, squares: s=15, w=100pm, TV=40000, stars: 
s=25, 10=130 /im, 7V=25000, triangles: s=20, w=10 pm, iV=15000. 
The dashed line is a numerical solution of equation ^6) for the 3D 
lattice using the density of states given by eq. 

For the general case, we use a low-temperature expansion 
similar to the Sommerfeld expansion for free electrons. We 
integrate eq. by parts and expand R(E) = J p(E)dE in 
a power series around E = p,, analogous to the method de- 
scribed in l2lll . We obtain the number of doubly occupied 
lattice sites to be 

P {E)dE + an(k B T) n n P _\ 1 | e =„ (8) 

-°° n=l ^ 

with the coefficients a n defined by 

f°° x n ( d 1 \ 

an= L^\-d-*wTw) dx (9) 

and given by a n = (2 — 2 2 ~ n )( t (n) for even values of n and 
o-n = — a-n-i f° r °dd values of n. The Sommerfeld expansion 
yields N = P (E)dE+Y,Zi MW 2 "^#U 
jH with the coefficients b 2n = ^ (-*-^) dx = 

a% n , which represent the even terms in the expansion l|8}. The 
analytic expression for the fraction of doubly occupied lattice 
sites is then given by 

n=0 




FIG. 3: Fraction of doubly occupied lattice sites 712 in a three- 
dimensional optical lattice as a function of the lattice depth and the 
temperature obtained by diagonalizing Hamiltonian Q. The particle 
number is N = 15000 per spin state and the waist of the lattice laser 
is w — 70 pm, according to the situation in |7J|. 

The dependence of the fraction of doubly occupied sites 
on temperature in three dimensions is stronger than in two di- 
mensions but in both cases we find for k B T <C Ep a predomi- 
nantly linear dependence with a slope of order unity. For com- 
parison we compute the fraction of doubly occupied lattice 
sites numerically. In figure^ we compare the results in two 
dimensions with the exact analytical prediction of equation 
0, which shows excellent agreement. In figure [2J? the nu- 
merical results for a three-dimensional optical lattice are com- 
pared with the expression dl 2I >. which is valid for kgT <C Ep. 
We also find very good agreement in this situation. These re- 
sults are also in agreement with the numerical results of 1 8]. 
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When the system is not in the low-tunnelling regime we 
leave the range of validity of the analytic calculation. To quan- 
tify this limit, we have numerically studied how the fraction 
of doubly occupied lattice sites varies with the potential depth 
of the lattice (Fig. [3}. We have chosen the experimentally 
realized values of N — 15000 per spin state and a waist of 
w = 70 /im [7] as parameters. From the calculations we ob- 
serve that the analytical solution, which is independent of the 
lattice depth, is good for a potential depth V > 12 E r . In 
this low-tunnelling regime, n 2 depends on temperature only. 
For smaller values of the lattice depth numerical calculations 
are necessary to relate a measured fraction of doubly occupied 
lattice sites to the temperature. 

Let us now compare our results with a recent experiment 
10. A noninteracting, two-component Fermi gas with an ini- 
tial temperature in the harmonic trap of T/Tp — 0.25 is 
loaded into a three-dimensional optical lattice. In a deep opti- 
cal lattice (V = 15 E r ) molecules are created by adiabatically 
sweeping a magnetic field across a Feshbach resonance and 
the molecule fraction is determined to be n 2 = 0.43. From 
our theory we calculate a temperature of T/Tp = 0.46 in the 
lattice. Provided that the sample in the lattice is in thermal 
equilibrium (which, however, has not been verified experi- 
mentally), this indicates a temperature increase of approxi- 
mately a factor two as compared to the initial situation. This 
is expected for adiabatic loading of the optical lattice due to 
entropy conservation in the loading process. 

Finally, we would like to address the weakly interacting 



two-component Fermi gas in an optical lattice. We take the 
on-site interaction between different spin states of strength U 
into account using a mean-field approach 1 18], which requires 
that kp\a\ <C 1. The occupation number v a {fi) per spin state 
a at a lattice site in the low-tunnelling regime is given by 
vM) = (expKVir^ + U^in^ir^fi-ri/kBTl + l)- 1 . 
This equation is solved self-consistently. The results show that 
the fraction of doubly occupied sites changes with the strength 

of U: For repulsive interaction ri2 = J2i V ]{ r i) v l{ r i)/N 
is diminished and for attractive interaction n 2 is increased. 
However, the interaction induced effect is smaller than the 
temperature effect. For kp\a\ = 0.05, w = 70 /im, N — 
15000 (the parameters of figure^ and a lattice depth of 25 E r 
we observe a change of double occupancy by approximately 
2%, which is small as compared to the effect of temperature. 

In conclusion, we have investigated a Fermi gas in an op- 
tical lattice with an underlying harmonic trapping potential. 
Our results have immediate consequences for ongoing experi- 
ments with ultracold fermions in optical lattices. We find that 
the Fermi gas is adiabatically heated by approximately a fac- 
tor two when the gas is loaded into the lattice. Moreover, we 
calculate how the occupation of an optical lattice depends on 
temperature and find that the fraction of doubly occupied lat- 
tice sites for a noninteracting, two-component Fermi gas is a 
sensitive quantity for thermometry in the lattice. 
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